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Preface 


The  goal  of  this  study  was  to  determine  design  guidelines  for  the  SP-100 
space  reactor  core  cooling  system  during  a  loss-of-coolant  accident  (LOCA).  These 
guidelines  were  determined  by  analysis  of  generic  fuel  assemblies  common  to  the 
SP-100  designs.  The  assemblies,  each  assembly  contains  36  individual  fuel  pins 
and  one  cooling  bayonet,  were  modeled  through  the  use  of  the  TSAP  computer 
program.  This  procedure  identified  regions  of  influence  for  cooling  within  each 
fuel  assembly.  Once  the  region  of  influence  is  identified  for  the  assemblies,  an 
entire  SP-100  reactor  core  can  be  analyzed  to  determine  the  maximum  cooling  sys¬ 
tem  necessary  to  ensure  survival  of  LOCA  conditions. 

Since  this  study  uses  generic  design  features  of  the  SP-100  reactor  core,  it 
may  be  extended  to  the  competing  SP-100  derivative  reactor  designs.  This  study 
is,  therefore,  an  analysis  and  design  tool  for  continued  studies  into  the  SP-100 
program. 

My  appreciation  and  thanks  go  to  my  faculty  advisor,  Maj  D.E.  Beller  ,  for 
his  guidance.  Also,  thanks  go  to  Maj  J.A.  Lupo  for  his  computer  expertise,  and  to 
Lt  Col  R.F.  Tuttle  for  his  support  and  enthusiasm.  Special  thanks  goes  to  the 
staff  of  Los  Alamos  National  Laboratories  for  sponsoring  my  work  and  providing 
assistance.  In  particular,  I  would  to  thank  C.  Bell,  F.  Biebl,  J.  Elson,  and  J. 
Boudreau  for  their  technical  guidance  and  without  whose  help  this  work  would 
not  have  been  completed. 

Finally,  I  would  like  to  express  ray  appreciation  to  my  family,  Linda,  Emily, 
and  Donna,  for  standing  by  me  and  supporting  me  for  the  last  eighteen  months. 
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Abstract 

The  purpose  of  this  study  was  to  determine  design  guidelines  for  the  SP-100 
space  reactor  core  cooling  system  after  a  loss-of-  cool  ant  accident.  SP-100  design 
parameters  were  obtained  from  Los  Alamos  National  Laboratory.  The  Thermal 
Systems  Analysis  Code  (TSAP)  calculated  the  temperatures  within  the  fuel  assem¬ 
blies  as  a  result  of  the  fuel  decay  heat.  TSAP  is  a  lumped-parameter  network 
analysis  code  capable  of  performing  radiative  and  conductive  heat  transfer 
analysis. 

The  reactor  core  was  assumed  to  void  of  coolant  instantaneously  following  a 
LOCA.  The  reactor  core  model  consisted  of  individual  fuel  pin  assemblies  con¬ 
taining  36  fuel  pins  surrounding  a  central  cooling  channel.  This  central  cooling 
channel,  or  bayonet,  is  a  secondary  cooling  loop  within  the  reactor  core.  The  bay¬ 
onet  cooling  is  a  safety  feature  designed  to  keep  the  core  from  reaching  tempera¬ 
tures  at  which  the  uranium  dissociates  from  the  uranium-nitride  fuel.  TSAP  cal¬ 
culated  the  fuel  pin  temperatures  due  to  decay  heat  and  identified  the  maximum 
temperatures  throughout  the  decay  heat  transient.  The  performance  of  the  bay¬ 
onets  within  a  generic  reactor  core  was  compared  to  an  actual  design.  Design  gui¬ 
dance  was  established  based  on  the  performance  of  the  bayonets  in  the  generic 
core. 
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AN  ANALYSIS  OF  HEAT  TRANSFER  AFTER  LOSS  OF  PRIMARY 
COOL.YNT  IN  THE  SP-100  REACTOR  SYSTEM 


I.  Introduction 


1.1  Background 

The  Department  of  Defense  (DOD),  Department  of  Energy  (DOE),  and  the 
National  Aeronautics  and  Space  Administration  (NASA)  are  jointly  developing 
space  reactor  technology  (1:232).  This  technology  development  is  in  response  to 
military  and  civilian  space  efforts  which  demand  the  use  of  small,  durable,  high- 
power  systems.  Examples  of  these  systems  include  the  Strategic  Defense  Initia¬ 
tive,  space-based  radar  and  navigational  systems,  and  the  orbiting  space  station. 
The  energy  requirements  of  these  and  future  efforts  in  space  require  greater  power 
and  a  longer  operational  lifetime  than  is  currently  available  from  conventional 
auxiliary  power  systems.  Currently,  the  program  chartered  by  the  DOD,  DOE, 
and  NASA  for  developing  space  nuclear  power  technology  is  the  SP-100  Program 
(1:232). 

Some  of  the  SP-100  Program  objectives  include  (1:232): 

1.  Ensure  aerospace  nuclear  safety. 

2.  Define  performance  limits  and  satisfy  military  and  civilian  applications  for 
space  nuclear  power. 

3.  Advance  technology  in  a  100-kilowatt  space  nuclear  power  system  to  support 
military  and  civilian  missions. 

4.  Evaluate  needs  for  multimegiwatt  space  nuclear  power  systems  and  perform 
long  range  technology  development  efforts. 
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5.  Effectively  transition  from  technology  assessment  and  advancement  to 
engineering  development  and  production  for  space  reactor  missions. 

The  safety  issue  is  one  of  the  biggest  problems  faced  by  the  SP-100  Program  (14). 
Three  phases  to  the  safety  issue  have  been  defined  by  DOE  (14).  These  phases 
are: 

1.  During  launch. 

2.  In  orbit,  during  operational  mode. 

3.  In  orbit  during  retirement. 

Los  Alamos  National  Laboratory  (LANL)  is  actively  involved  in  assessing  the 
safety  of  each  phase.  At  this  time,  LANL  believes  the  most  critical  safety  problem 
is  designing  the  reactor  to  withstand  a  loss-of- cool  ant  accident  while  it  is  in  orbit 
(14).  General  Electric  completed  a  study  (8)  showing  the  vulnerability  of  the  SP- 
100  to  impact  from  meteorites  and  debris  left  in  orbit  from  other  missions. 

Impact  by  debris  on  the  reactor  vessel  can  rupture  the  vessel  wall  causing  a  break 
in  the  primary  core  cooling  system.  A  loss-of-coolaut  accident  (LOCA)  can  occur 
when  the  primary  coolant  piping  breaks  and  the  coolant  can’t  circulate  through 
the  core.  Even  though  safety  rods  reduce  the  reactor  power,  there  is  a  continuing 
supply  of  heat  from  decaying  fission  products  in  the  core.  This  condition  can 
cause  the  fuel  within  the  core  to  increase  in  temperature  until  it  begins  to  melt 
and  damages  the  structural  properties  of  the  fuel  rods.  Once  the  structural 
integrity  of  the  reactor  is  damaged,  retiring  the  system  to  a  high  orbit  may 
become  difficult.  There  is  also  a  high  probability,  after  the  threshold  temperature 
is  exceeded,  that  molecular  uranium  will  be  released  in  a  low  earth  orbit  as  the 
uranium  dissociates  from  the  fuel.  Therefore,  before  deployment,  a  complete  ther¬ 
mal  analysis  must  be  performed  to  ensure  adequate  core  cooling  after  a  LOCA. 


1.2  Problem 


This  study  will  investigate  the  hypothesis  that  a  generic  SP-100  reactor  core 
can  be  used  as  a  design  and  analysis  tool  for  other  SP-100  designs.  Figure  1.1  is  a 
picture  of  a  typical  SP-100  reactor  core.  It  is  made  up  of  many  fuel  assemblies 
and  internal  structural  material.  The  internal  structure  gives  the  core  its  strength 
by  holding  the  assemblies  in  place  and  adding  a  conductive  heat  transfer  path. 

.Figure  1.2.  shows  .a.  generic  fuej. assembly  totally  void  of  internal  structure,  except 
for  the  fuel  cladding. 

If  the  fuel  is  allowed  to  melt  the  internal  structure  will  begin  to  be  attacked 
by  the  molten  fuel  and  chemical  reaction  from  free  molecular  uranium.  The 
uranium  dissociates  from  the  uranium-nitride  fuel  at  about  2000  K.  By  keeping 
the  fuel  temperature  below  2000  K  the  core  will  not  melt  and  the  internal  struc¬ 
ture  will  remain  intact.  The  problem  row  becomes  one  of  keeping  the  core  below 
the  threshold  temperature. 

To  understand  how  to  keep  the  core  below  the  threshold,  this  study  will 
investigate  the  behavior  of  the  individual  fuel  assemblies  which  make  up  the  core. 
The  assemblies  will  be  void  of  internal  structure  and  will  differ  from  one  another 
only  by  the  fuel  rod  linear  power.  By  placing  these  assemblies  in  a  core 
configuration,  as  shown  in  figure  1.1,  the  problem  of  core  cooling  after  a  LOCA 
can  be  analyzed. 

1.3  Scope 

This  study  will  be  concerned  with  the  SP-100,  100-kilowatt  reactor  and  its 
derivatives.  Specific  design  features  considered  are:  initial  fuel  pin  temperatures, 
coolant  temperatures,  emissivity  of  the  fuel  pins,  and  the  fuel  linear  power.  The 
reactor  parameters  are  given  in  Appendix  A.  The  Thermal  Systems  Analysis  Pro¬ 
gram  (TSAP)  (19)  will  be  used  along  with  the  decay  heat  data  in  Appendix  E  to 
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find  the  maximum  fuel  temperature.  From  the  maximum  fuel  pin  temperatures 
design  and  analysis  guidelines  can  be  devised  by  considering  the  regions  below  the 
2000  K  threshold,  the  temperature  at  which  the  uranium  dissociates  from  the  fuel 
(10)- 

The  following  design  parameters  will  be  used  along  with  the  data  in  Appen¬ 
dix  A  to  model  the  reactor  core: 

■  pin  diameter 

■  emissivity  of  fuel  pins 

■  coolant  temperature 

■  initial  fuel  pin  temperature 

■  fuel  pin  linear  power 

1.4  Assumptions 

Assumptions  used  in  this  study  include: 

1.  Only  a  1-centimeter  horizontal  slice  of  the  core  mid-plane  is  modeled. 

2.  Fuel  pins  are  characterized  by  the  initial  pin  temperature,  pin  thermal  capa¬ 
citance,  and  the  pin  decay  heat  generation  time  history. 

3.  Heat  pipes  and  bayonet  cooling  rods  are  modeled  by  constant  temperature 
only. 

4.  Heat  transfer  between  fuel  pins  is  radiative. 

5.  Heat  transfer  is  in  the  radial  direction  and  negligible  in  the  axial  direction. 

6.  Loss  of  coolant  is  instantaneous,  with  respect  to  the  time  scale  of  the  tran¬ 
sient  coolant  loss,  and  completely  voids  the  core  of  coolant . 

Details  concerning  these  assumptions  can  be  found  in  the  following  chapter. 


=  0.838  cm 
=  0.05  -  0.5 
=-  800-  1500  K 
=  1000  -  1400  K 
=  30  -  100  watts/cm 


1.5  Approach 

The  approach  of  this  study  is  to  develop  a  generic  model,  void  of  internal 
structure,  of  the  SP-100  reactor  core.  The  model  will  consist  of  an  individual, 
generic  fuel  assembly  void  of  conductive  heat  transfer  paths.  This  assembly  will 
consist  of  a  central  bayonet  surrounded  by  hexagonal  rows  of  fuel  pins.  From  the 
design  parameters,  changes  to  the  model  will  be  made  to  determine  the  efL:t  on 
the  pin  temperatures.  Data  will  be  taken  on  the  number  of  fuel  pin  rows  which 
remain  below  the  2000  K  threshold  temperature.  This  procedure  will  continue 
until  all  the  design  parameters  have  been  exercised  through  the  ranges  of  interest 
defined  by  LANL. 

After  accumulating  data,  the  trends  indicated  by  the  results  will  be  applied 
to  an  SP-100  reactor  core  design.  Design  guidance  procedures  will  be  developed  as 
the  reactor  is  analyzed  using  the  results  from  the  generic  model.  In  addition, 
recommendations  for  future  work  will  be  included  in  the  final  chapter. 


II.  Theory 


2.1  Introduction 

The  hypothesis  that  a  generic  SP-100  reactor  core  can  be  used  to  determine 
design  guidelines  for  a  complete  SP-100  reactor  core  must  be  based  on  basic 
scientific  and  engineering  principles.  This  chapter  demonstrates  those  principles 
used  in  the  analysis  of  the  generic  reactor  core  and  allows  the  application  of  the 
data  discussed  in  chapter  3  to  any  SP-100  configuration  with  the  parameters  exhi¬ 
bited  in  Appendix  A. 

2.2  Computer  Code 

A  lumped-parameter  computer  code  called  Thermal  Systems  Analysis  Pro¬ 
gram  (TSAP)  was  used  for  the  calculations.  It  is  a  thermal  network  code  which 
treats  the  heat  transfer  as  an  electrical  circuit  with  the  fuel  pins  and  bayonets 
representing  homogeneous  nodes  (11).  Between  any  two  black  surfaces  1,2  the 
radiant  heat  exchange  rate  is  given  by 

01,2  “  A  1F12  1^61  “  Eb2  j  (2-1) 

This  equation  is  an  analog  to  Ohm’s  law,  where  the  quantity  of  transfer  q; 

the  potential  driving  force  Ekx—Ek2\  and  the  thermal  resistance,  - - - ;  have 

a\F\,2 

electrical  conterparts  /,  AP,  and  R  respectively  (2:251).  The  emissive  power  of  a 
gray  body  is  given  by  the  Stefan-Boltzmann  Law  (2:437): 

Ek  -  £  a  T4  (2.2) 

where  e  represents  the  emissivity  of  the  surface.  Emissivity  is  a  material  property 
and  can  be  determined  experimentally.  Figure  2.1  shows  the  typical  model  ele¬ 
ments  incorporated  in  TSAP,  and  the  electrical  analog  heat  transfer  paths. 


2.1 


(a)  Fuel  pins  and  assembly  walls  are  characterized  by  Initial  temperatures, 
thermal  capacitance,  and  a  heat-generation  profile* 

(b)  Type  1  and  Type  2  conductors  are  radiative  couplers. 

(c)  Type  3  conductors  are  conductive. 


(d)  Type  4  conductor*  are  radiative  except  for  a  variation  In  Configuration  4 
where  they  are  conductive. 


(e)  A  heat  pipe  Is  characterized  by  a  fixed  temperature  and  a 
pipe  by  a  thermal  capacitance  and  an  Initial  temperature. 


Figure  2.1  TSAP  Model  Elements  (Ref  11) 
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2.3  Computer  Model 


A  computer  model  of  a  SP-100  reactor  core  was  developed  as  an  input  file  to 
TSAP.  The  model  was  void  of  any  internal  structure,  which  may  have  contri¬ 
buted  a  conductive  path  in  the  model;  the  core  contained  only  fuel  pins  and  bay¬ 
onet  type  cooling  rods.  In  order  to  simplify  the  analysis,  the  outer  boundary  of 
the  model  was  taken  as  adiabatic.  This  results  in  a  50100  K  higher  maximum 
temperature  than  the  case  of  allowing  the  heat  to  escape  into  space.  However,  if 
the  reactor  can  be  cooled  at  this  slightly  elevated  temperature,  it  will  be  within 
the  safety  threshold  of  2000  K  in  an  operational  mode.  In  addition,  the  model 
used  (Appendix  B)  takes  advantage  of  symmetry,  so  then  only  one  small  segment 
of  the  reactor  core  need  be  modeled.  The  actual  configuration  shown  in  figure  1.2 
can  be  seen  to  be  a  subassembly  of  the  fuel  pin  assemblies  shown  in  figure  1.1. 

Appendix  B  contains  all  the  information  necessary  for  TSAP  to  calculate  the 
maximum  fuel  pin  temperatures  following  a  LOCA.  The  "NODE  DATA"  section 
contains  the  data  on  each  individual  fuel  pin  or  bayonet  node.  Bayonet  nodes  are 
assumed  to  be  the  same  size  as  fuel  pins  and  are  actually  a  generic  cooling  surface. 
The  bayonet  could  be  a  heat  pipe,  a  cold-finger,  or  any  single  point  cooling 
mechanism.  All  operational  cooling  elements  are  characterized  by  the  first  line 
contained  under  the  line  NODE  DATA.  The  negative  sign  indicates  a  constant 
temperature  node,  the  number  101  is  simply  a  bookkeeping  notation  for  the  node 
number.  Bayonet  temperature,  heat  capacity,  and  linear  power  complete  the 
input  for  the  bayonet.  Since  the  bayonet  is  a  cooling  device,  the  heat  capacity 
and  power  level  is  zero.  However,  if  the  bayonet  has  failed  then  the  model  may  be 
changed  by  removing  the  negative  sign  to  indicate  a  node  of  changing  tempera¬ 
ture,  in  this  case  the  bayonet  heat  capacity  must  be  supplied  for  TSAP  to  calcu¬ 
late  the  temperature  change  within  the  bayonet. 

Fuel  pin  data,  the  lines  following  the  bayonet  data,  contains  the  pin  number, 


the  pin  initial  temperature  in  degrees  kelvin,  the  fuel  pin  thermal  capacitance  in 
joules  per  degree  kelvin,  and  the  fuel  pin  linear  power  in  watts  per  centimeter  of 
pin  length.  Thermal  capacitance  is  discussed  in  section  2.4  of  this  study.  Since 
the  fuel  pin  data  is  not  preceded  by  a  negative  sign,  the  fuel  pin  temperature  is 
not  held  constant,  and  is  therefore  varied  according  to  the  heat  transfer  paths 
open  to  the  fuel  pin.  The  heat  transfer  path  data  is  contained  under  the  "COND 
DATA"  section  of  the  TSAP  input  file  contained  in  Appendix  B. 

The  heat  transfer  paths  or  conduction  data  can  be  identified  by  TS.AP  as 
either  radiative  or  conductive.  If  a  conduction  data  line  begins  with  a  negative 
sign  then  TS.AP  used  a  radiative  heat  transfer  path,  if  not  then  a  conductive  path 
is  used.  The  first  number  specifies  the  conductor  identification  number,  the  next 
two  numbers  refer  to  the  two  nodes  the  conductor  connects.  These  node  numbers 
refer  to  the  bayonets  or  fuel  pins  from  the  "NODE  DATA"  section  of  the  input. 
The  last  entry  in  "COND  DATA"  is  the  conductance  parameter  given  in  watts  per 
degree  kelvin  to  the  fourth  power.  This  parameter  is  discussed  in  the  next  sec¬ 
tion. 

The  individual  nodes  and  their  conductance  paths  have  been  defined,  now  the 
heat  sourace  must  be  included.  For  the  case  of  a  LOCA  the  only  heat  source  is 
the  decay  heat  of  the  fission-products.  The  "ARRAY  DATA"  section  includes  the 
time  history  of  the  decay  heat,  and  the  percent  of  total  thermal  power  the  decay 
heat  represents.  The  first  line  is  the  time  entry  in  seconds,  and  the  second  line 
contains  the  decay  heat  data. 

Other  sections  of  the  TSAP  input  file  include  the  "CNTL  DATA"  which 
specifies  time  steps  for  TSAP  to  use  in  calculations,  and  the  remainder  of  the 
input  describes  the  calculation  of  the  heat  flux  for  use  within  the  TSAP  calcula¬ 
tions. 


2.4  TSAP  Pin-to-Pin  Heat  Transfer  Modeling 


The  TSAP  code  requires  a  thermal  capacitance  input  for  each  node  and  a 
thermal  conductance  input  for  each  conductive  path  connecting  nodes  (11).  For 
the  thermal  capacitance,  C,  we  have: 

C  =  C~m  (2 

where 


m  =  p  V 


and 


P  =  £/,  P. 

i 


= 


£/,  Pi  c, 

i 


c,  is  the  specific  heat  of  material  i 


C„  is  the  average  specific  heat  of  node 


F,  is  the  volume  fraction  of  material  i 


m  is  the  node  mass 


V  is  the  node  volume 


p,  is  the  density  of  material  i 
p  is  the  average  density  of  the  node 


Therefore,  for  a  fuel  pin  node  the  thermal  capacitance  is 

CP'n  =  T  r*ad  h  (  (PCP  \lad  +  (PCP  Lrr+  (PCP  )/„e/  ) 
For  the  bayonet  or  heat,  pipe  nodes 

Cb  =  Cpt  pb  7r  h  \rcjad  -  rj ) 

TSAP  uses  the  conductance,  G ,  between  nodes  t  and  j  such  that 

1  -  c  (r,  -  t,  ) 


for  conductive  heat  transfer,  and 


=  g  (r,4  -  r/) 


for  radiative  heat  transfer,  where 


q  is  the  heat  transfer  between  nodes  t  and  j 


T,  is  the  temperature  of  node  j 


T}  is  the  temperature  of  node  j 


The  following  is  a  derivation  of  the  conductance  parameter  G  between  pins  i 
and  j  (9).  Each  pin  is  assumed  to  have  identical  emissivities,  thermal  conductivi¬ 
ties,  densities  and  specific  heats.  The  fuel  pins  are  cylindrical  and  are  positioned 
in  an  array  on  a  triangelar  pitch.  Additionally,  the  pins  are  assumed  infinite  in 
length;  each  pin  is  assumed  to  radiate  only  to  its  six  nearest  neighbors;  the  inter¬ 
stitial  gaps  between  pins  are  small  compared  to  some  unit  fuel  pin  length,  and 
each  pin  is  represented  by  a  single  node  and  a  single  temperature  (9).  This 
hypothesises  that  temperature  gradients  within  a  given  fuel  pin  are  small  com¬ 
pared  to  radiative  temperature  drops  between  pins. 


Figure  2.2  is  used  in  this  derivation  to  define  the  surfaces  in  question.  From 
the  figure  the  heat  transfer  between  »  and  j  is  defined  as  twice  the  heat  from  sur¬ 
face  1  to  surface  3. 

=2?i-3  (2-7) 

Assume  pin  i  and  pin  j  have  equal  temperatures,  then  the  heat  transfer 
between  surface  1  and  surface  3  is: 

<?  1-3  =  ~2  Q 1-2  3  (2-8) 

Combining  equations  2.7  and  2.8  yields 

?i->  “  <7  1-2,3  (2-9) 


Looking  at  the  expression  for  two  gray  bodies  that  only  see  each  other; 


a,<t 


where 


(2.10) 


A  j  is  the  area  of  surface  1 
A  2  3  is  the  area  of  conbined  surfaces  2  and  3 

F i_2  3  is  the  view  factor  from  surface  1  to  the  combined  surfaces  2  and  3 

T,  is  the  temperature  of  pin  i 

Tj  is  the  temperature  of  pin  j 

e  is  the  material  emissivity 

rr  is  the  Stefan-Boltzmann  constant 


2.7 


Combining  equations  2.6,  2.9,  and  2.10 


A  2.3  =  2  1 


F  1-23  —  1.0 

For  a  fuel  pin  of  diameter  D  and  length  L 


Simplifying 


G 


—  7T  D  L  <7 

3 


e 


This  expression  was  derived  using  the  assumption  of  equal  pin  temperatures 
for  adjacent  pins.  Elson  in  reference  9  derives  this  equation  as  a  3-body  radiant 
exchange  problem. 


2.5  Axial  Conduction 


Throughout  the  study,  heat  transfer  in  the  axial  direction  has  been  assumed 
negligable  compared  to  the  radiation  heat  transfer  in  the  radial  direction.  To 
demonstrate  that  this  assumption  is  valid  for  reactors  of  this  size,  the  axial  heat 
transfer  and  radial  heat  transfer  are  compared. 


In  the  axial  direction  the  heat  flux  rate  per  rod  is  found  from 


(2.1  T 


Where  K  is  the  thermal  conductivity  of  the  material.  The  radial  heat  flux  rate 
per  rod  is  determined  by  equation  2.10.  Using  infromation  from  Appendix  A  and 
reference  18,  the  axial  conductive  heat  transfer  rate  per  rod  is  0.14  watts  and  the 
radial  radiative  heat  transfer  rate  per  rod  is  27.0  watts.  These  calculations  were 
based  on  the  smallest  temperature  difference  considered  for  the  SP-100  design,  200 
K.  This  temperature  difference  minimized  the  T4  effect  of  the  radiative  heat 
transfer.  When  larger  changes  in  temperature  are  allowed  (300-  SO  K)  the  radia¬ 
tive  heat  transfer  rate  grows  faster  than  the  conductive  heat  transfer  rate. 

For  this  case  the  ratio  of  axial  conduction  versus  radial  radiation  is  on  the 
order  of  5xlO-3.  Therefore,  neglecting  the  effect  of  axial  conductive  heat  transfer 
by  the  fuel  rods  is  justified. 


III.  Results  of  the  Parametric  Investigation 


3.1  TSAP  Procedures 


The  parameters  of  emissivity,  fuel  pin  temperature,  coolant  temperature,  and 
fuel  pin  linear  power  were  varied  over  the  ranges  listed  in  Chapter  I.  These  values 
were  input  into  the  TSAP  code  as  discussed  in  Chapter  II  and  in  reference  II.  A 
listing  of  the  sample  input  is  in  Appendix  B. 

TSAP  uses  the  node  and  conduction  data  to  create  a  thermal  map  of  the 
reactor  core  region.  The  decay  heat  time  history  in  Appendix  E  is  used  to  calcu¬ 
late  the  temperature  of  each  node  as  the  decay  heat  continues  to  drop  off  in  time. 
A  new  thermal  map  of  the  core  region  is  generated  for  each  time  step  indicated  in 
the  input  file  under  the  "CNTRL  DATA"  entries.  TSAP  will  also  generate  a  tem¬ 
perature  history  of  the  core  region  as  a  function  of  time,  identifying  the  hottest 
node  within  the  region  and  the  time  at  which  this  temperature  occurs.  From  this 
updated  temperature  data,  the  number  of  fuel  pins  remaining  below  the  2000  K 
threshold  are  identified  and  the  number  of  hexagonal  rows  containing  these  pins 
can  be  determined.  In  addition,  the  total  amount  of  pin  linear  power  cooled  by  a 
bayonet,  located  in  the  central  part  of  the  mapped  core  region,  can  be  calculated. 

Two  parameters  which  are  important  in  a  LOCA  are  the  fuel  pin  emissivity 
and  the  coolant  temperature  in  the  secondary  cooling  system.  The  data  displayed 
on  the  following  graphs  are  arranged  to  show  cooling  trends  as  one  is  held  con¬ 
stant  and  the  other  is  varied.  Another  important  parameter  in  this  analysis  is  the 
Region  of  Influence.  The  features  of  this  parameter  are  the  actual  number  of  hex¬ 
agonal  rows  contained  in  the  region  and  the  power  (watts/cm)  which  these  rows 
represent.  Figures  3.1  -  3.  4  demonstrate  these  features. 


Figure  3.1  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temperature  =  1000  K) 
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Figure  3.2  Region  of  Influence  (watts/cm)  versus  Linear  Power 
(fuel  Temperature  =  1000  K) 
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Figure  3.3  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temperature  =  1000  K) 
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Figure  3.4  Region  of  Influence  (watts/cm)  versus  Linear  Power 
(fuel  temperature  =  1000  K) 
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3.2  Results 


Figures  3.1  and  3.2  show  the  Region  of  Influence  for  the  case  of  a  fixed  emis¬ 
sivity,  e  —  0.3.  Figures  3.3  and  3.4  are  fixed  at  a  coolant  temperature  of  800  K. 

All  of  these  figures  are  for  a  reactor  fuel  assembly  with  a  fuel  pin  temperature  of 
1000  K.  The  linear  fuel  pin  power  is  on  the  x-axis  and  the  Region  of  Influence  is 
displayed  on  the  y-axis.  Symbols  for  coolant  temperature  and  emissivity  on  the 
top  graph  on  each  page  also  hold  for  the  other  graph  on  the  page. 

For  a  fixed  emissivity,  the  region  of  influence,  in  number  of  rows  cooled,  docs 
not  change  much  (1  row)  with  coolant  temperature  (700  K)  at  the  same  linear 
power  level.  However,  if  the  coolant  is  fixed  (800  K)  and  the  emissivity  is  allowed 
to  vary,  the  number  of  rows  cooled  varies  from  4  to  6  at  30  watts/cm  to  no  varia¬ 
tion  at  the  100  watts/cm  power  region.  This  demonstrates  that  if  the  cladding 
emissivity  changes  an  additional  heat  removal  mechanism  can  be  gained  or  lost 
depending  on  the  change  in  the  emissivity.  If  the  emissivity  increases  with  the 
fuel  pin  age,  due  to  roughening  effects,  an  older  core  has  a  better  chance  of  surviv¬ 
ing  a  LOCA  than  a  fresh  core. 

From  Figure  3.1  at  a  fuel  pin  linear  power  of  30  watts/cm,  for  an  emissivity 
of  0.5,  a  coolant  at  800  K  will  cool  6  hexagonal  rows  of  fuel  pins.  The  assumption 
here  is  that  the  coolant  is  located  at  the  center  of  the  hexagonal  fuel  array.  At  a 
linear  power  of  100  watts/cm,  the  same  conditions  can  only  cool  2  hexagonal 
rows.  For  systems  considering  high  power  densities,  the  emissivity  of  the  clad  is 
an  important  parameter  to  maximize  for  inherent  safety  during  a  LOCA. 

Figure  3.4  shows  the  power  contained  in  the  region  of  influence  held  under 
the  threshold  temperature.  The  apparent  discontinuities  are  due  to  the  discrete 
nature  of  the  definition  of  region  of  influence.  The  region  of  influence  assumes  a 
hexagonal  heat  transfer  geometry  due  to  fuel  assembly  manufacturing  and  design 
techniques  (15).  However,  the  heat  transfer  path  is  not  hexagonal  in  nature. 


3.4 


TSAP  produces  a  map  of  the  temperatures  within  a  given  core  and  the  results, 
not  part  of  this  study,  indicate  a  different  geometric  pattern.  The  sloping  pattern 
in  figure  3.4  indicates  that  the  coolant  is  capable  of  cooling  at  3780  watts/cm  but 
at  100  watts/cm  per  pin  the  same  conditions  cool  only  1800  watts/cm  (for 
e  «■  0.5).  This  indicates  that  the  fuel  pins  in  the  hexagonal  array  at  100  watts/cm 
are  probably  overcooled. 

The  remaining  graphs  are  included  in  Appendix  F.  Additional  data  was 
taken  for  an  emissivity  of  e  =  0.05.  At  this  low  emissivity  the  coolant,  at  any 
temperature  in  the  region  800  -  1500  K,  could  cool  at  the  most  1  row  at  power 
levels  between  30  -  50  watts/cm.  Above  this  power  level  all  the  rods  melted 
regardless  of  the  coolant  temperature.  If  the  emissivity  drops  for  any  reason  (eg., 
lithium  plating),  the  core  will  not  survive  a  LOCA  under  the  current  designs. 

This  data  is  included  in  the  tables  in  Appendix  C. 
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IV.  Discussion  and  Application 


4.1  Discussion 

The  graphs  in  the  previous  chapter  show  a  strong  dependency  on  the  fuel  pin 
emissivity.  In  the  cases  where  the  initial  fuel  pin  temperature  and  the  coolant 
temperature  are  held  constant,  only  the  emissivity  is  allowed  to  vary  for  a  given 
power  level.  Regions  of  influence  for  maintaining  the  core  below  the  threshold 
value  vary  considerably  for  the  different  power  levels.  The  change  in  the  number 
of  rows  cooled,  from  2  in  the  low  power  regions  (30  -  50  watts/cm)  to  no  variation 
in  the  high  power  region  (  >  50  watts/cm),  indicates  that  high  power  density 
cores  have  an  inherent  problem  with  LOCA.  The  lower  coolant  temperatures  in 
the  bayonets  does  not  relieve  the  problem. 

The  graphs  of  region  of  influence  in  terms  of  linear  power  cooled  show  a  wide 
variation.  This  indicates  that  the  cooling  bayonets  can  cool  more  fuel  pins  in  a 
region;  however,  the  region  is  not  hexagonal  in  geometry.  This  leads  to  the 
assumption  that  the  cooling  bayonets  are  overcooling  the  fuel  pin  assemblies. 

If  a  reactor  is  designed  with  the  cooling  bayonets  placed  in  the  center  of  the 
fuel  assemblies,  then  the  number  of  hexagonal  rows  expected  to  remain  below  the 
threshold  is  given  by  the  Regions  of  Influence  (Rows).  The  Region  of  Influence 
(watts/cm)  indicate  how  much  linear  power  the  bayonet  can  support  before  the 
threshold  is  reached.  This  parameter  is  a  constant  for  a  given  bayonet.  However, 
due  to  the  definition  of  the  Region  of  Influence  and  the  method  for  allowing  row 
survival  the  bayonet  is  capable  of  cooling  a  larger  region.  This  region  is  not  hex¬ 
agonal  in  geometry. 

The  method  of  allowing  row  survival  is  based  on  the  criterion  that  all  pins  in 
the  row  must  remain  below  the  2000  K  threshold  to  be  counted.  This  is  based  on 
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the  reactor  core  design  criteria  and  core  fabrication  techniques.  (Fuel  assemblies 
are  usually  placed  in  a  triangular  pitch  for  neutronic  and  heat  transfer  purposes 
(15)).  Since  a  triangular  pitch  is  used  in  fabrication  the  fuel  assemblies  form  hex¬ 
agonal  shapes. 

4.2  Application 

To  demonstrate  the  design  guidance  from  the  graphs  a  reactor  core  was 
chosen  from  one  of  the  preliminary  SP-100  design  candidates.  A  3-safetv-rod 
configuration  was  chosen  to  demonstrate  the  applicability  of  the  preceding  data. 
The  3-safety-rod  design  (Appendix  D)  failed  at  a  temperature  of  2507  K.  The 
design  features  a  cooling  rod  temperature  of  800  K,  fuel  pin  temperature  of  1320 
K,  and  an  emissivitv  of  0.3.  The  power  varied  from  51  watts  /cm  at  the  outside 
edge  to  63  watts/cm  at  the  central  part  of  the  core.  The  coolant  is  allowed  to 
flow  around  the  safety  rod  and  is  maintained  at  a  constant  temperature. 

The  figures  in  Appendix  F  contain  data  for  a  reactor  with  slightly  cooler  fuel 
pins  (1250  K)  but  the  same  emissivity  and  coolant  temperature.  From  this  data, 
it  can  be  shown  that  for  the  given  power  levels  from  the  graphs,  the  largest  region 
of  influence  of  the  coolant  is  2  hexagonal  rows  of  pins.  As  a  result,  it  can  be 
deduced  that  this  reactor  is  undercooled  after  LOCA  and  will  fail.  However,  if  a 
bayonet  is  placed  at  the  center  of  each  assembly,  the  required  cooling  is  achieved. 
Each  fuel  assembly  contains  a  centra!  fuel  pin  surrounded  by  36  other  pins  in  a 
hexagonal  array.  This  establishes  a  region  of  influence  of  3  for  each  central  pin  if 
that  pin  were  replaced  with  a  bayonet.  The  reason  the  reactor  achieved  a  tem¬ 
perature  below  threshold  is  due  to  the  conductive  paths  now  available  and  the 
overcooling  of  the  pin  within  the  2  row  region  of  the  bayonet. 

When  the  3-safety-rod  model  was  run,  using  the  cooling  system  suggested  by 
figure  3.9,  the  threshold  temperature  was  never  achieved  (maximum  was  1879  K). 
In  addition,  the  conductive  paths  within  the  3  rod  model  transferred  the  heat  to 


4.2 


the  outside  of  the  core  faster  and  the  peak  core  temperature  occured  much  later 
then  the  decay  heat  peak. 

From  this  application  two  important  facts  become  apparent.  First,  the 
figures  in  chapter  3  can  be  used  to  quickly  gage  the  appropriate  number  of  cooling 
bayonets  required  in  each  region.  Second,  conductive  paths  help  the  cooling  pro¬ 
cess  by  driving  the  peak  core  temperature  to  a  later  time  than  the  decay  heat 
peak  temperatures.  Since  application  of  the  data  from  figure  3.0  resulted  in  the 
core  maintaining  a  temperature  of  1879  K  instead  of  2507  K,  the  concept  of  using 
the  generic  model  is  valid.  The  conductive  paths,  which  are  not  included  in  the 
generic  modei,  act  to  enhance  the  cooling. 

Further  studies  need  to  be  performed  on  the  fuel  pin  materials  to  determine 
the  effects  surface  roughness  and  age  have  on  the  emissive  properties.  If  some 
coolant  remains  in  the  pores  of  the  cladding,  the  emissivity  of  the  fuel  cladding 
will  decrease.  However,  if  the  fuel  cladding  becomes  rougher  with  age  the  emis¬ 
sivity  will  increase. 
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V.  Summary  and  Recommendations 


5.1  Summary 

This  study  performed  a  thermal  analysis  of  a  SP-100  reactor  core,  void  of  all 
internal  structure,  after  a  loss-of- cool  ant  accident.  The  purpose  was  to  determine 
the  maximum  temperatures  achieved  in  a  region  and  to  adjust  the  region  so  the 
maximum  temperature  did  not  exceed  2000  K,  the  dissociation  temperature  of  the 
uranium  from  the  uranium-nitride  fuel. 

Results  of  this  study  produced  a  set  of  design  graphs  which  can  be  applied  to 
determine  the  initial  number  of  bayonets  needed  in  the  cooling  system.  When 
applied  to  a  core  which  initially  failed  after  LOCA  due  to  inadequate  cooling;  the 
new  cooling  system  (dictated  by  the  data)  enabled  the  reactor  core  to  maintain 
temperatures  well  below  the  threshold. 

The  results  also  indicate  a  strong  dependence  on  the  fuel  pin  emissivity.  This 
dependence  has  a  stronger  influence  on  the  reactor  cooling  system  than  the  tem¬ 
perature  of  the  coolant.  The  data  in  Appendix  C  and  Appendix  F  show  that  the 
coolant  temperature  may  vary  as  much  as  500  K  before  a  difference  is  noted  in  the 
region  of  influence.  However,  if  the  emissivity  changes  by  the  same  percentage  a 
larger  variation  can  be  seen,  especially  at  low  power  densities. 

In  addition  to  the  material  studies,  further  work  needs  to  be  completed  to  fill 
in  the  graphs  in  this  study.  As  mentioned  earlier,  due  to  the  definition  of  region 
of  influence,  the  exact  behavior  between  the  data  points  cannot  be  easily 
identified.  For  example,  is  the  region  of  influence  a  step  function  of  the  power,  or 
is  there  an  empirical  relationship  somewhere  in  the  data?  The  answer  to  this 
question  will  greatly  enhance  the  utility  of  the  data  in  this  study. 


Overall,  the  generic  design  study  showed  that  a  successful  design  can  be 
accomplished  using  the  data  obtained.  The  data  in  this  study  will  result  in  over¬ 
cooling  a  specific  reactor  core  design  due  to  the  nature  of  the  Region  of  Influence 
(Rows).  Emissivity  also  plays  a  major  role  in  the  cooling  after  a  LOCA.  The 
exact  behavior  of  the  emissivity  in  the  fuel  pins  can  be  used  to  gage  the  number  of 
bayonets  needed  to  remain  below  the  threshold.  Aging  and  lithium  plating  are 
contributing  factors  to  the  emissivity  changes  within  the  core  and  have  a  direct 
bearing  on  the  emissive  properties  of  the  pins. 

5.2  Recommendations 

As  mentioned  earlier,  the  emissivity  of  the  cladding  is  critical  to  the  heat 
transfer  of  the  core  during  a  LOCA.  Further  studies  need  to  be  completed  on  the 
emissivity  changes  due  to  plating  and  aging.  Surface  roughness  tests  should  be 
performed  to  simulate  these  effects.  The  geometric  pattern  of  preferential  heat 
transfer  within  the  fuel  assemblies  should  be  studied  to  determine  the  extent  of 
overcooling  offered  by  the  data  presented  in  this  study.  Since  the  hexagonal 
region  of  influence  results  in  overcooling  the  core,  a  better  method  of  analyzing 
the  heat  transfer  in  the  core  may  be  found  within  the  natural  heat  transfer 
geometry. 
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APPENDIX  A:  SP-100  Design  Data  For  TSAP  Analysis 


reactor  thermal  po war  (W) 

m 

2.50000a+06 

number  of  fuel  pin  assemblies 

m 

3.10000*^1 

number  of  pins  par  assembly 

m 

3.70000s*01 

eora  height  (a) 

m 

3.00380a-01 

peek-to-average  power  factor 

m 

1.44G00e*O0 

fual  pallat  outar  radius  (a) 

m 

3.20000a-03 

fual/linar  gap  outar  radius  (a) 

m 

3.27500a-03 

liter  outar  radius  (a) 

m 

3.53000a-03 

elsd  outar  radius  (a) 

m 

4.19000a-03 

fual  dansity  (kg/m3) 

m 

1.35140e«-04 

llnar  density  (kg/m3) 

m 

2.10000e+04 

clad/structura  dansity  (kg/a3) 

m 

8.40000e+03 

fual  spacifle  hast  (J/kgX) 

m 

2 . 52000e+02 

llnar  specific  hast  (J/kgX) 

m 

1.54000e+02 

clad/structura  spacifle  hast  (J/kgX) 

m 

2.75000»H)2 

duct  outar  flat-to-flat  distanca  (a) 

m 

5.74500a-02 

duct  wall  threknasa  (a) 

m 

1.24500S-03 

duct  wall  conductivity  (*/mX) 

m 

6.56800a+01 

hast  pipa  wall  thickness  (a) 

- 

7.00000a-04 

axial  no da  height  (m) 

- 

1.00000--02 

clad/structura  aaissivity 

- 

3.00000a-01 

fual  pin  capacitance  (J/X) 

m 

1.64162a-00 

heat  pipa  capacltanea  (J/K) 

m 

3.90140e-01 

duct  half-wall  capacltanea  (J/X) 

m 

4.«6622e-0l 

paak  pin  powar  (¥) 

- 

6.27248**01 

vi 


I 

: 

t 

i 

i 

i 

* 

I 

» 

<• 

. 

■> 

■. 
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pla-to-pin  conduct ancs  (W/K4)  • 

adga  pin-to-duct  conductanct  (X/X4)  • 

comas  pin-to-duct  conductanca  (W/K4)  - 

half  duct-to-half  duct  conductance  (rid-W/X4)» 
half  duct-to-half  duct  conductanca  (eond-x/K)  - 


5.52828e-13 

9.85977a-13 

7.2l244a-13 

1.65933e-12 

4.93064a-02 


Figure  A.l  SP-100  Design  Information  (Courtesy  LANL) 
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APPENDIX  B:  TSAP  Input  File 


AFIT  GENERIC  REACTOR  CORE  MODEL  (Ti=1000,Tcf=800,P=10,E=.3) 
NODE  DATA 
-101,800.0,0.0,0.0 
201,1000.0,0.8208,5.0 
301,1000.0,0.8208,5.0 
302,1000.0,0.8208,5.0 
401,1000.0,1.6415,10.0 
402,1000.0,0.8208,5.0 
501,1000.0,0.8208,5.0 
502,1000.0,1.6415,10.0 
503,1000.0,0.8208.5.0 
601,1000.0,1.6415,10.0 
602,1000.0,1.6415,10.0 
603,1000.0,0.8208,5.0 
701,1000.0,0.8208,5.0 
702,1000.0,1.6415,10.0 
703,1000.0.1.6415,10.0 
704,1000.0,0.8208,5.0 
801,1000.0,1.6415,10.0 
802,1000.0,1.6415,10.0 
803,1000.0,1.6415,10.0 
804,1000.0,0.8208.5.0 
001,1000.0,0.8208,5.0 
902,1000.0,1.6415,10.0 
903,1000.0,1.6415,10.0 
904,1000.0,1.6415,10.0 
905,1000.0,0.8408,5.0 
1001,1000.0,1.6415,10.0 
1002,1000.0,1.6415,10.0 
1003,1000.0,1.6415,10.0 
1004,1000.0,1.6415,10.0 
1005,1000.0,0.8408,5.0  ENDD  COND  DATA 
-5001, 101, 201, 2.7642E-13 
-5002,201 ,301 .5.5283E-13 
-5003,201 ,302,2.7642E-13 
-5004,301,302,5.5283E-13 
-5005,301 ,401 ,5. 5283E- 1 3 
-5006,302,401,5.5283E-13 
-5007,302,402,2.7642E-13 
-5008,401 ,402,5.5283E-13 
-5009, 401, 501, 5.5283E-13 
-5010,401,502,5.5283E-13 
-5011,402, 502, 5.5283E-13 
-5012, 402, 503, 2.7642E-13 
-5013,501 ,502,5.5283E-13 
-5014, 502, 503, 5.5283E-13 
-5015, 501, 601, 5.5283E-13 
-5016, 502, 601, 5.5283E-13 
-50 1 7 , 502 ,602 .5.5283E- 1 3 


v 
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i: 


-5018,503,602,5.5283E-13 

-5019, 503, 603, 2.7642E-13 

-5020,601,602,5.5283E-13 

-5021, 602, 603, 5.5283E-13 

-5022,601,701 ,5.5283E-13 

-5023,601 ,702,5.5283E-13 

-5024 ,602,702,5.5283E- 1 3 

-5025, 602, 703, 5.5283E-13 

-5026, 603, 703, 5.5283E-13 

-5027, 603, 704, 2.7642E-13 

-5028, 701, 702, 5.5283E-13 

-5029, 702, 703, 5.5283E-13 

-5030,703,704 ,5.5283E- 1 3 

-5031, 701, 801 ,5.5283E-13 

-5032, 702, 801, 5.5283E-13 

-5033, 702, 802, 5.5283E-13 

-5034, 703, 802, 5.5283E-13 

-5035, 703, 803, 5.5283E-13 

-5036, 704, 803, 5.5283E-13 

-5037, 704, 804, 2.7642E-13 

-5038,801, 802, 5.5283E- 13 

-5039, 802, 803, 5.5283E-13 

-5040,803,804 .5.5283E- 1 3 

-5041,801,901 ,5.5283E-13 

-5042,801, 902, 5.5283E-13 

-5043, 802, 902, 5.5283E-13 

-5044,802, 903, 5.5283E-1 3 

-5045,803,903,5.5283E-13 

-5046, 803, 904, 5.5283E-13 

-5047,804 ,904 .5.5283E-1 3 

-5048, 804, 905, 2.7642E-13 

-5049,901 ,902,5.5283E-13 

-5050,902,903,5.5283E-13 

-5051, 603, 904, 5.5283E-13 

-5052, 904, 905, 5.5283E-13 

-5053, 901, 1001, 5.5283E-13 

-5054, 902, 1001, 5.5283E-13 

-5055,902, 1002, 5.52S3E  13 

-5056, 903, 1002, 5.5283E-13 

-5057, 903, 1003, 5.5283E-13 

-5058, 904, 1003, 5.5283E-13 

-5059, 904, 1004, 5.5283E-13 

-5060,905, 1004 ,5.5283E- 1 3 

-5061,905, 1005, 2.7642E-13 

-5062, 1001, 1002, 5.5283E-13 

-5063, 1002, 1003, 5.5283E-13 

-  5064, 1003,1 004, 5.5283E- 13 

-5065, 1004, 1005, 5.5283E-13  ENDD  CNTL  DATA 

IPRO,0, 

NL  OP, 5000, 

DCNV,0.02, 

ISTA,1, 

TSTP.IO., 

TSTA,0., 


TEND,  14401., 

•  TIME,0., 

IVAB,1»  ENDD  ARRY  DATA  1 
0.0, 60., 600., 3600., 36000., 8.64E+05,8.64E+06,END  2 
.0742, .0401, .0247, .0138, .00658, .00256, .000867, END  ENDD  VABl  DATA 
C 

PER=TERPL1(A1  ,A2, TIME, TIME,  l.,l.) 

Q201  =  5.000  *  PER 
Q301  =  5.000  *  PER 
Q302  =  5.000  *  PER 
Q401  =  10.000  *  PER 
Q402  =  5.000  *  PER 
Q501  =  5.000  *  PER 
Q502  —  10.000  *  PER 
Q503  =  5.000  *  PER 
Q601  =  10.000  *  PER 
Q602  =  10.000  *  PER 
Q803  =  5.000  *  PER 
Q701  =  5.000  *  PER 
Q702  =  10.000  *  PER 
Q703  =  10.000  *  PER 
Q704  =  5.000  ♦  PER 
Q801  =  10.000  *  PER 
Q802  =  10.000  *  PER 
Q803  =  10.000  ♦  PER 
Q804  =  5.000  *  PER 
Q901  =  5.000  *  PER 
Q902  =  10.000  *  PER 
Q903  =  10.000  *  PER 
Q904  =  10.000  *  PER 
Q905  =  5.000  *  PER 
Q1001  =  10.000  *  PER 
Q1002  =  10.000  *  PER 
Q1003  =  10.000  *  PER 
Q1004  =  10.000  *  PER 
Q1005  =  5.000  *  PER  ENDD  VAB2  DATA 
C 

EF(TIME  .LT.  7201.)THEN 

ITIME— NINT(TIME) 

IF(ITIME  .GT.  3601)THEN 
NPRT=1 

ELSEIF(ITIME  .GT.  3509)THEN 
TS  TP =3600. 

ELSEIF(ITIME  .GT.  1801)THEN 
NPRT=3 

ELSEIF(ITIME  .GT.  601)THEN 
NPRT=2 

ELSEIF(ITIME  .GT.  599)THEN 
TSTP=600. 

ELSEIF(ITIME  .GT.  421)THEN 
NPRT=3 

ELSEIF(IT£ME  .GT.  301)THEN 
NPRT=2 

ELSEIF{ITIME  .GT.  61)THEN 
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NPRT=2 

ELSEIF(ITIME  .GT.  59)THEN 
TSTP=60. 

ELSE 

TSTP=10. 

NPRT=3 
END  IF 
END  IF 

CALL  OUTPUT(NDTOTL,Q,N,TIME,EROR,LOOP.4) 

CALL  QFLOW(NDTOTL,NCTOTL,T,G,IZl,IZ2,IZ3,M,IPRO)  ENDD 
VAB3  DATA 

C  ENDD  USER  DATA 

C******USe  this  statement  in  vab2  to  get  output  of  heat  flux  at  end 
C******  of  each  time  step 
C* 

C* 

c* 

C*  CALL  QFLOW(NDTOTL,NCTOTL,T,G,IZl,IZ2,IZ3,M,IPRO) 

C* 

c* 

c* 

SUBROUTINE  QFLOW(ND,NC,T,G,IZl ,IZ2,IZ3,M,IPRC) 


DIMENSION  T(_l ),G(1  ),IZl ( 1  ),IZ2(1  ),IZ3(  1  ),M(  1 ) 
DIMENSION  QFLO( 30000) 

C 

TABS=  0.0 

IF  (IPRO  .NE.  1)  TABS—  0.0 

IACUM=0 

DO  100  1=1, ND 

C -  number  of  conductors  connected  to  node  I 

IZ=1Z1(1) 

IF  (IZ  .EQ.  0)  GO  TO  100 
IF  (IZ  .LT.  0)  1Z=-IZ 
DO  200  J=1,IZ 
JI=IACUM+J 

C -  what  node  it  connects  to 

IZP=IZ2(JI) 

C — conductor  number  -  TSAP 
IZG=IZ3(JI) 

C -  radiation 

IF  (IZG  .LT.  0)  THEN 
IZG=-IZG 

C —  actual  conductor  number 
TI=T(I)+TABS 
T1ZP=T  IZPJ+TABS 
QFLO(M(izg))=G(IZG)*(TIZP**4-Tl“4) 
ELSE 

QFLO(M(izg))=G(IZG)*(T(IZP)-T(I)) 

END  IF 


*4 


r  ,type  *,.,?/Mde’,,*l/.?,nnec^  toMZP/with  conductor  #\m(izg) 

£  type  *,’t(i)=',T(i),  T(izp)=’,T(izp),’  g=»\GflZGV 

$low ( ’'n°80 » )•  )  =\Qflo(nog(ji$ 


;  »  vw.—  (izp),  g: 

200  CONTilS^OW(’,n0g(ji)’  }  -^MnogOijj 

iacum  =  iacum  •+•  IZ 
100  CONTINUE 
WRITE(6,1) 

iformat(//,4X/***  HEAT  FLUX  ***’) 

2^ORM.\T(2X,6(lX,l5,2X,Ell.4)) 

RETURN 
END  ENDD  ENDD 
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APPENDIX  Cj  Table  of  Results  for  the  100  kWe  SP-100  from  TSAP 
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2 

.30 
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60 
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3 

.50 

6 
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3780 

4 

50 

.05 

1 

6 
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5 

.30 

3 

36 
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6 

.50 

4 

60 
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7 

70 

.05 

0 

0 

0 
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18 
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APPEND  DC  E:  SP-100  Fiaaion-Product  Decay  Heat 


Time  (Sec 


%  Full  Power 


7.42 

6.82 


6.49 


5.95 


10.0 

20.0 


60.0 


5.45 


4.90 


4.01 


120.0 


3.47 


300.0 


2.86 


600.0 


2.47 


1200.0 


1.91 


10800.0 


1.38 


36000.0 


0.963 


86400.0 


0.499 


259200.0 


0.366 


864000.0 


0.256 


2592000.0 


0.164 


8640000.0 


0.0857 


31536000.0 


0.0311 


94608000.0 


0.0128 


315360000.0 


0.00744 


E.l 


APPENDIX  F:  Graphical  Display  of  Data 


The  data  presented  here  is  the  graphical  form  of  the  data  in  Appendix  C 
except  for  the  exclusion  of  the  data  corresponding  to  an  emissivity  of  0.06.  This 
data  was  discussed  earlier  and  contains  no  useful  design  information  other  than  to 
avoid  that  emissivity.  The  other  curves  however,  are  displayed  and  should  prove 
useful  to  reactor  designers  investigating  LOCA  conditions. 

The  symbols  represented  on  the  top  graph  on  each  page  also  correspond  to 
the  symbols  on  the  bottom  of  the  same  page. 
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Figure  F.l  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.2  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.3  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.4  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.5  Region  of  Influence  (Rows]  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.6  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.7  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.8  Region  of  Influence  (W/CM)  versus  Linear  Power 

(fuel  temp.  =  1000  K) 
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Figure  F.9  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.10  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.ll  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F. 12  Region  of  Influence  (VV/CM)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.13  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.14  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  —  1000  K) 
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Figure  F.15  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.16  Region  of  Influence  (VV/CNl)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 


F.e 


FIXED  COOLANT  TEMPERATURE  (1300  k) 


Linear  Tower  (wnlls/om) 


Figure  F.17  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.18  Region  cf  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1000  K) 
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Figure  F.10  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  letup.  =  1000  K) 
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i'  igure  F  -20  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  —  1000  K) 
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Figure  F.21  Region  of  Influence  (Rows)  versus  Linear  Pow'er 
(fuel  temp.  =  1250  K) 


FIXED  COOLANT  TEMPERATURE  (000  k) 


6 

^  5 

5  * 

e  " 

» 

—  3 
8  " 
£  S 

3 


a 

& 

1 


l  - -  -1 


1  *  t 


o 


□ 


O 


u 


m 


-j  _ i - j - j  . 


J  _  .  — i  _ l 


H  II  l.l. 

I.lnrnr  Power  (wnlls/ein) 


-'o  Mu 


Figure  F.22  Region  of  Influence  (VV/CM)  versus  Linear  Power 
(fuel  tcrnp.  =  1250  K) 
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Figure  F.23  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.24  Region  of  Influence  (W/CM)  versus  Linear  Pow'er  ; 

(fuel  temp.  =  1250  K)  | 
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Figure  F.25  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.27  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.28  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1250  K) 
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Figure  F.29  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.30  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.31  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.32  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.33  Region  of  Influence  (Rows)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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Figure  F.34  Region  of  Influence  (W/CM)  versus  Linear  Power 
(fuel  temp.  =  1400  K) 
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